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Abstract: We present preliminary results of Particle-In-Cell simulations of magnetic turbulence produc- 
tion by isotropic cosmic-ray ions streaming upstream of supernova remnant shocks. The studies aim at 
testing the MHD predictions by Bell [1, 2] of a strong amplification of short- wavelength nonresonant 
wave modes and at studying the subsequent evolution of the magnetic turbulence and its backreaction 
on cosmic ray trajectories. The detailed knowledge of the upstream turbulence properties is crucial to 
ascertain all aspects of the shock acceleration process - the transport properties of cosmic rays, the shock 
structure, thermal particle injection and heating processes. An amplification of magnetic field would also 
facilitate the acceleration of particles beyond the "knee" in the cosmic-ray spectrum. Our kinetic ap- 
proach is particularly suited to address the backreaction on the cosmic rays, and it allows us to test Bell's 
prediction of the eventual formation of extended filamentary structure in the cosmic-ray distribution and 
also to arrive at a reliable estimate of the total saturation magnetic-field level. The parameters chosen for 
the simulations are favorable for the rapid excitation of purely growing modes. We confirm the generation 
of the turbulent magnetic field due to the drift of cosmic-ray ions in the upstream plasma, but show that 
the growth rate of the field perturbations is much slower than estimated using the MHD approach and 
the amplitude of the turbulence saturates at about SB/ B ~ 1. The magnetic field also remains below 
equipartition with the upstream plasma. 



Introduction 

Supernova remnant shock waves are prime candi- 
dates for the acceleration sites of Galactic cosmic 
rays. However, no firm evidence has been found 
so far for hadron production in SNRs, which chal- 
lenges our understanding of these sources. Par- 
ticle confinement near the shock is supported by 
self-generated magnetic turbulence and thus the 
detailed knowledge of the properties of the turbu- 
lence is crucial to ascertain all aspects of the ac- 
celeration process - transport properties of cosmic 
rays, the shock structure, thermal particle injection 
and heating processes. Of particular interest is the 
question of how efficiently and with what proper- 
ties electromagnetic turbulence is produced by en- 
ergetic particles at some distance upstream of the 
forward shocks. As the flow convects that turbu- 
lence toward and through the shock, the turbulent 
magnetic field structure at and behind the shock is 
shaped by the plasma interactions ahead of it. 



We investigate the properties of magnetic turbu- 
lence upstream of the shocks of young SNRs. Cos- 
mic rays accelerated at these shocks stream rela- 
tive to the upstream plasma with the shock veloc- 
ity, with the highest energy particles propagating 
the farthest from the shock and forming a nearly 
monoenergetic quasi-isotropic population at some 
distance upstream. Such a system has been re- 
cently studied by Bell [1] with a quasilinear MHD 
approach. Bell noted that the current carried by 
drifting cosmic rays excites non-resonant, nearly 
purely growing modes more efficiently than reso- 
nant Alfven waves, which he verified with MHD 
simulations. Later studies suggested that regions of 
low density and low magnetic field strength would 
be formed, thus resembling a filamentation insta- 
bility [2]. Here we use Particle-In-Cell simula- 
tions to verify the existence of the non-resonant 
modes in the presence of kinetic particle effects 
and explore its properties such as nonlinear evo- 
lution, saturation, and particle heating. 
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Figure 1: Temporal evolution of the total energy density in fields (a) and particles (Ucr, Ui, Ug) and the 
magnetic field Um (b) normalized to the initial total kinetic energy in the system Uko- Time is in units of 
the growth rate ■jmax for the most unstable mode with wavelength Xmax- B\ and E\ indicate the magnetic 
and electric field energy densities in the components perpendicular to the cosmic-ray ion drift direction. 
The total kinetic energy in cosmic-ray ions Ucr is rescaled by a factor of 10. 



Particle-in-cell experiment and results 

We use a modified version of the relativistic 3D 
code TRISTAN. In the simulations, an isotropic 
population of relativistic, monoenergetic cosmic- 
ray ions with Lorentz factor ^cr = 2 and num- 
ber density Ncr drifts with Vgh = 0.3c relative 
to the upstream electron-ion plasma and along the 
homogeneous magnetic field i?||o, carrying a cur- 
rent density jcR = eNcRVsh- The ions of the 
upstream medium have a thermal distribution with 
number density Ni, in thermal equilibrium with 
the electrons. The electron population with density 
Ng = Ni + Ncr contains the excess electrons to 
provide charge-neutrality and also drifts with Vd = 
VshNcR/Ne with respect to the upstream ions pro- 
viding a return current jret = —eN^Vd to balance 
the current in cosmic rays. The density ratio as- 
sumed in this study is NcR/Ni=l/3. The simula- 
tions have been performed using a 992x304x304 
computational grid with periodic boundary con- 
ditions and with a total of 1.5 billion particles. 
Cosmic-ray ions drift in the — a:-direction. The 
electron skindepth A^e = c/cOpe = 7 A, where 
ojpe is the electron plasma frequency and A is the 
grid cell size, and the ion-electron mass ratio as- 



sumed mi/nie = 10. The above parameter com- 
bination was chosen to resemble physical condi- 
tions in young SNRs and, according to quasiUn- 
ear calculations [1], is favorable for the rapid ex- 
citation of purely growing, short-wavelength (as 
compared with cosmic-ray ion gyroradii rcRg) 
wave modes. In the unstable wavevector regime 
1 < k^^rcRg < CVgf^/v\, the dispersion relation 
for the purely growing mode 

-vlkl + Cvlh/rcRg^O (1) 

gives the maximum growth rate jmax at wavenum- 
ber k^^rnax = C,vlj^l2v\rcRg, where k\\ is wavevec- 
tor parallel to Byo, va is the Alfven velocity, 
and C = NcRCjcR/NiVsh- Our parameters give 
C,'v1h/v\ ~ 826 and the most unstable mode with 
wavelength Xmax = 50A is in the middle of the 
unstable wavevector range. 

Figure 1 shows the temporal evolution of the mag- 
netic and electric field and particle kinetic energy 
densities. An initial (up to t ^ Inun-) f'lst growth 
of the turbulent fields is caused by a Buneman-type 
beam instability between slowly drifting electrons 
and the ambient ions, which leads to plasma heat- 
ing until the electron thermal velocities become 
comparable to the electron drift speed Vd. In real 
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Figure 2: Power spectra of the magnetic field component By in the direction parallel (a) and perpendicular 
(b) to the cosmic -ray ions drift direction at different times as given in units of the growth rate -fmax- The 
spectra were calculated in the 2D slice in the xz-plane at y/A — 200. 



SNRs the density ratio is likely much larger than 
the sound Mach number, and therefore the electron 
drift speed is smaller than the thermal velocity. The 
initial instability is thus solely a consequence of the 
initial conditions in our simulations, in particular 
the initial temperature of the background medium. 

The wave modes due to the instability caused by 
the cosmic -ray ions streaming in the upstream 
plasma start to emerge att^ '^Imax ^i^^ '^he mag- 
netic field turbulence is excited mainly in the com- 
ponents transverse to the cosmic-ray drift direc- 
tion. The dominant wave mode is oblique at an 
angle of about 75° to the s-direction (Fig. 2), and 
has the wavelength component « 50A, which 
is numerically similar to Xmax- The character of 
the mode is thus different than that predicted from 
the linear analysis by Bell [1, 2], which assumed 
that the most rapid growth occurs for wavenum- 
bers parallel to i?||o (Eq. 1). As seen in Fig. 2, 
the dominant scale in X± grows as the turbulence 
develops, but the parallel scales (A|| = A^;) remain 
roughly constant, so after about 207^^^ the struc- 
tures have similar size parallel and perpendicular 
to the drift. Comparison with Fig. 1 reveals that 
at those late times 5B± is of the order of i3||o, so 
the system becomes magnetically nearly isotropic. 
The drift of the cosmic -ray particles is still kinet- 
ically dominant, which suggests an important role 



of the homogeneous magnetic field component for 
the evolution of the magnetic turbulence. 

The spatial structures of the turbulent magnetic 
field and accompanying ambient plasma distribu- 
tion are shown in Fig. 3 for t « 107^^3.. The 
perpendicular field is concentrated around regions 
of low background plasma density. The cosmic- 
ray distribution is homogeneous, and therefore 
the magnetic field is associated with a strong net 
cosmic -ray current flowing inside ambient-plasma 
voids, and the field lines circle around the center 
of the cavities according to Ampere's law. The 
formation of the cavities in the plasma is due to 
the jret X B force, which accelerates the upstream 
plasma away from the center of the voids, and 
also causes the cavities to expand. This is visi- 
ble as the evolution towards larger scales shown in 
Fig. 2b. The growth of the cavities and the subse- 
quent merging of the adjacent plasma voids lead to 
a compression of the plasma between cavities and 
also to an amplification of the magnetic field. Be- 
cause the magnetic field has a preferred orientation 
around each cavity, the magnetic field lines cancel 
each other in the space between the voids (Fig. 3). 

Compared with Bell's MHD studies, the growth of 
the magnetic perturbations in our kinetic model- 
ing is much slower The initial growth rate {t ^ 
(2 — 6)^max) of the perpendicular field turbulence 
is only 1/5 of ''jmax, and becomes smaller during 
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Figure 3: The magnitude and direction of the perpendicular magnetic field component — {By+Bfy/^ 
(a) and the ambient ion density Ni (b) in the plane perpendicular to the cosmic-ray ion drift direction at 
x/A — 500 and t « IQj^l^^. B± is normalized to the amplitude of the homogeneous field B\^q, and Ni to 
the initial ambient ions density. The electron distribution follows that of ambient ions. 



the later evolution. Moreover, the amplitude of the 
turbulent component never considerably exceeds 
the amplitude of the regular field and the field 
growth saturates when SB±/B\\q « 1. Note, that 
although for t > 257"^^. the field structures be- 
come larger than the size of the simulation box, the 
magnetic field saturation level is not influenced by 
the boundary conditions, which we have confirmed 
using larger-box 2D simulations, that well repro- 
duce the results of the 3D experiment. The field 
growth is also accompanied by heating of the am- 
bient plasma, and the field remains below equipar- 
tition with the upstream medium. The cosmic -ray 
density distribution remains roughly unchanged till 
the end of the simulation. 

Summary and conclusions 

The process of the production of magnetic field 
turbulence by cosmic-ray ions drifting upstream of 
SNR shocks has been studied here using PIC sim- 
ulations. Turbulent field is indeed generated in this 
process, but the growth of magnetic turbulence is 
much slower than estimated using the MHD ap- 
proach, and the amplitude of the field perturbations 



saturates at approximately the amplitude of the ho- 
mogeneous upstream field. The energy density 
in the turbulent field is also always much smaller 
than the plasma kinetic energy density. This sug- 
gests that the efficiency of magnetic field gener- 
ation in this mechanism may not be sufficient to 
account for the high magnetic fields observed in 
some young SNRs and also leaves open the ques- 
tion of the ability of diffusive particle acceleration 
at SNR shocks to produce particles with energies 
beyond the "knee" in the cosmic -ray spectrum. 

The simulations have been performed on Columbia at 
NAS and Tungsten at NCSA which is supported by the 
NSF. The work was supported by MNiSW during 2005- 
2008 as research project 1 P03D 003 29. 
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